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THE TRANSPORT OF VORTICITY THROUGH FLUIDS 
IN TURBULENT MOTION* 
(In the light of the Prandtl and Taylor theories) 

By C. Ferrari 



SUMMARY 



The author studies the problem of the transport of 
vorticity or of momentum in the light of the Taylor and 
Prandtl theories which he "briefly reviews. He shows how 
the formulas of Prandtl could he brought into agreement 
with experimental results in those cases where they agree 
with the principle of statistic similitude of Karman, and 
particularly in the problem of the distribution of veloc- 
ity and temperature in the wake of a heated cylindrical 
obstacle. He shows that when the formulas are extended to 
two-dimensional motion with streamlines whose curvature is 
not zero, they lead to unsatisfactory results and that in 
this case the formulas differ completely from those de- 
rived from the principle of similitude when the latter is 
applied either to the configuration of disturbed flow or 
to the distribution functions of the turbulent velocities. 
He then examines the relations of this problem of trans- 
port of motion with the theory of Mattioli, which appears 
susceptible of some observations. After pointing out the 
difficulty of obtaining a satisfactory theory of turbu- 
lence based 'on the concept of transport and deduced by the 
methods of classical mechanics, he indicates the reasons 
therefor and shows finally how the problem may find a so- 
lution by applying the methods of statistical mechanics 
according to the theory of G-ebelein. 

1» The essential defining characteristics of the tur- 
bulent motion of a fluid is the well-known irregular fluc- 
tuat ions , both in magnitude and d irect ion , of the velo c ity 
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at each point, corresponding to an energetic mixing of the 
mass whereby a fluid particle may, in the course of its 
motion, occupy any position whatever in the field. The 
fluctuations/Of the" modulus and argument of the velocity 
vector about their mean value, though "both irregular, are 
not, however, independent. There exists between them a 
statistical correlation, in consequence of which the mean 
time value of the product of two velocity components along 
any two perpendicular directions is different from zero. 
There result upon each element of surface immersed in the 
fluid virtual stresses perfectly analogous to the viscous 
stresses produced by thermal molecular agitation. To char- 
acterise this turbulent agitation of the mass Prandtl (ref- 
erence 1) and Taylor, (reference 2) have independently in- 
troduced the concept of mixing length or "Mischungsweg" , 
denoted by I, which is quite analogous to the molecular 
mean free path X considered in the kinetic gas theory, 
and which may therefore be defined as the path normal to 
the lines of flow which the particles can trace out and 
still maintain their individuality; that is, without as- 
suming the physical characteristics of the medium in which 
they are immersed. But whereas Prandtl and his collabor- 
ators assume that throughout the path I each fluid parti- 
cle maintains its momentum, Taylor objects that the in- 
stantaneous differences of pressure may cause the velocity 
of the displaced particle to vary and that therefore it is 
more logical to assume that it is the vorticity, upon which 
the instantaneous local pressure distribution has no ef- 
fect, that is maintained constant. 

2. In order to understand these two concepts more . 
clearly, let us limit our considerations to. two-dimensional 
fields of motion in which the flow lines are exactly or 
approximately straight lines parallel to the X axis, along 
which the characteristics of the motion may be assumed con- 
stant. If u' and v' are the two components of the ve- 
locity due to the. turbulent agitation, the transport of 
the momentum across an element dx is equal to 

T xy dx = - p u'v' dx (l) 

where u'v' denotes the mean time value of the product of 
u» and v 1 ; the increment of momentum per unit volume com- 
municated to the fluid layer of thickness dy is there- 
fore 
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and therefore if the mean, notion takes place in the direc- 
tion of the X axis with a pressure gradient 3p/9x, the 
equation of motion is " " ' " " • ■ ■ ■ • - • • • •■• 

- I E = - J" (3) 
p 3x - dy 

pn the other hand, considering the motion from the 
Lagrahgian point of view, the fluctuations in velocity- 
appear as a consequence of the transfer of fluid particles 
from layer to. layer, so that if l» is the distance, at a 
given instant, "between the layer from which the particle 
comes to that, which it occupies, and. if. it is assumed with 
Prandtl that within the distance V the velocity remains 
constant, then 

u< = - l» J2 (4) 
dy 

where U denotes the velocity of the motion. From (3) 
there results 

p 3x dy V V dy J KO) 

Equation (5) has further "been transformed "by Prandtl 
"by making the plausible assumption that the correlation 
"between u' and v' and therefore "between I 1 and v 1 
is constant over the whole field, so that it is possible 
to put 



l 1 v' 

in which c is a constant and l" the mean square varia- 
tion of l' , or, according to the definition given above, 
the mean "mixing length." Considering next two fluid par- 
ticles which, moving from the layers of height y .+ I 1 and 
y - V i meet at the layer of height y; they will approach 
each other or move away from each other with a relative 
velocity 2u' and will therefore induce in the fluid a 
velocity v' which, on account of the continuity of the 
fluid mass itself, should : be of the same order of magni- 
tude as u' : that is, . v'. '= ail' ' and therefore, 

S7t* =;o x J^ r 0l i |a ( 6 ) 



= c Jv 2 Jv** = c I Vv» 



where c x is likewise a constant. 
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By means of (6) equation (5) i s' transformed into 



in which 



I d£ = A. J 
p dx dy 1 



dU 
dy 



is included in I. 



dU 
dy 



(?) 



3. To deduce the increment of momentum communicated 
per unit volume' to the element of thickness dy under the 
hypothesis that not the momentum but the vorticity remains 
constant during the transport of the fluid particles, 
Taylor considers the equation of motion of a perfect fluid 
which, under the assumption that the mean motion is uni- 
form and parallel to the direction of the X axis, assumes 
the expression 



p 3x 



(8) 



in which 



2 W a*; 



is the vorticity; and therefore, in the fields of motion 
considered above the increment of motion due to the tur- 
bulent agitation is, according to Taylor 



3t 



xy 

ay 



(9) 



if and v' are the instantaneous values of 

cillations of the vorticity and of the component 
the velocity. There is then obtained the equati 
sponding to (3) 



the os- 
v' of 
quation corre- 



P dx 



(30 



in which the correlation between d)' and v' is brought 
about by the same causes that produce that between u' 
and v 1 in the theory of Prandtl. On the other hand, if 
l' has the meaning defined above, then according to Taylor, 



0) = 



1 

2 



d 3 U 



and therefore, 
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1 $p d 2 U 

- |~ = v i« j-s 

p ox dy 3 



(5-) 



or also, from (6) and the assumption of constant correla- 
tion between v 1 and I 1 



l Sp 
p 



= I 2 



dU 

dy 



d 2 U 



(?') 



4. Without entering for the present into the merit 
or justification of the criteria which led to equations 
(5) and (5«), these may simply be discussed, as has been 
indicated by Taylor (reference 3) and Page (reference 4) , 
on the basis of a comparison of the results to which they 
lead with those deduced from experiment. Now it is known 
that in the uniform flow about a flat plate, when I is 
set equal to cy as required from simple considerations 
of the homogeneity of the formulas and of 
tude , equation (7) leads immediately to a 
the velocity at the surface of contact of 
which fact has been brilliantly confirmed 

sical turbulence theory of Zarman (reference 5) and by the 
experimental investigations of Nikuradse (reference 6). 
There is, in fact, obtained by double integration of (7) 



dynamic simili- 
distribution of 
the plate itself 
both by the clas 



U = a log y + b (8) 
and from the tests of Nikuradse, 

= 5.5 + 5.75 log 7 -^H 

y t 0 / p u 



in which T Q i s the tangential stress at the plate. 

Equation (7 1 ), on the contrary, in the case which we 
are examining, does not give any significant result; in 

fact, if — — = 0, then either I should be equal to zero 

which case corresponds to nonturbul ent. mot ion , or f-— = 0 
j2« dy 

or, finally, Si_U- = o, and therefore the velocity at the 
dy 3 

surface of contact of the plate should vary either linear- 
ly or parabolically . 

Although the example just discussed appears to bear 
out the theory of Prandtl , a contrary result is obtained 
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if there is considered the phenomenon of turbulent diffu- 
sion in the wake "behind a cylindrical obstacle with axis 
of symmetry in the direction of the X axis. This problem 
has been studied theoretically and experimentally by 
Schlichting (reference 7), who has found that the velocity 
U along the X axis may be expressed by means of the for- 
mula 



u o - TJ _ x _ x/a f 



U 0 ' 

where U 0 is the velocity of the undisturbed stream. 

The velocity Y along y perpendicular to x, by 
the equation of continuity, is given by 



I_ = _ L f , r 



while the mixing length I is proportional to x l/2 , or 

I =ax l/2 . Uow if the wake is narrow and therefore y is 
small compared to x, then according to (7) there is ob- 
tained the equation of motion: 

Uo an = A \i^r an] = 2a2x an £ji = A x m | 2 u (9) 
0 ax ay L ay J ay ay 3 3 y ^ 



while according to (7») there is obtained: 



an = p 

dx 6 



v7 a!n = a s x an a!n = 4 _ an aln 



ay' 



ay ay- 



ay ay 



(9i) 



Equations (9) and (9 1 ) are formally identical, but 
the coefficient of turbulent transport resulting from (9) 
is double that contained in (9'); it follows that the dia- 
grams of velocity deduced from (9) and from (9') are iden- 
tical, but the results will not be identical if the trans- 
port coefficient calculated from (9) or (9') is applied to 
other problems intimately connected with this one. Let -as 
consider, for example, the temperature in the wake of the 
same obstacle which is assumed to be heated. The equation 
of the heat transport is 



U, 



a* 
ax 



_a_ 

By 



(TV 



ay j 



do) 
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and by comparing (10) with (9) it is immediately recog- 
nized that the latter is satisfied if we put 

T = b U (11) 

that is, the law of distribution of temperature in the 
wake of an obstacle should, according to Prandtl' s theory, 
coincide with that for the velocity. If, however, (11) is 
applied to (10), the latter is not transformed into (9») 
and therefore according to the theory of Taylor, the tem- 
perature diagram does not coincide with that for the ve- 
locity. Tor the velocity there is obtained 

U 0 - U «3/2, 2 



U 0 " u c 



= (1 - t*") 



where U c is the value of IT on the axis of the wake 

whose width is 2Y and % is equal to y/7 , for the 
t emperatur es 

m £.3/2 

i~ = 1 - £ 

*c 



where T c denotes the temperature on the axis. 

How the tests of A. Fage and of Falkner (reference 8) , 
carried out on two cylindrical heated obstacles of circular 
section and lenticular section, respectively, have shown 
an excellent agreement between the velocity and tempera- 
ture distribution agreeing with the theory of Taylor, and 
they particularly well bring out the difference between 
the temperature and velocity diagrams, respectively (fig. 
l). In this connection it should, however, be remarked 
that the similarity between the distribution diagrams of 
the temperature and velocity affirmed by the theory of 
Prandtl, has been experimentally confirmed by F. Flias 
(reference 9) for the flow about a flat plate. It there- 
fore appears that while the theory of transport of momen- 
tum is confirmed by experiments in problems of bound tur- 
bulence (at the contact of the solid wall) , at least for 
the very simple cases considered above, the theory of 
transport of vorticity gives better results in the prob- 
lems of free turbulence. Taylor and Prandtl attribute 
this singular behavior to the fact that in the case of 
bound turbulence, as demonstrated by. the experiments of 
Fage and Townend (reference 10), the perturbations which 
make the particles move from layer to layer are produced 
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essentially "by vortices arranged so that their axis is 
parallel to the wall of the obstacle and along the undis- 
turbed stream, so that the phenomenon of turbulent agita- 
tion is three-dimensional, whereas in free turbulence the 
transport of the fluid particles is due principally to 
vortices arranged with their axis normal to the flow line 
of the mean motion and to the plane of the motion. In the 
first case migrations of the particles can take place with 
out having the velocity influenced by the differences in 

local pressure; at any rate, Taylor observes that__the 

transport of vorticity is now given not only by 2w T z v' , 
but by 2(W^ v' - lO 1 ^ w 1 ) , where z is the direction nor- 
mal to the plane of the mean motion and w' and fo'y the 

component of the oscillation velocity and of the rotation, 
respectively, along the Z and Y axes. From the general 
laws of vortex motion of Helmholtz (see Lamb, Hydrodynam- 
ics, 4th ed., p. 197), Taylor deduces that the transport 
of vorticity in the case of three-dimensional perturba- 

3T ~. 3z dw z 

tions is given by v' — — u) z + I v< — - — , instead of by 

3c 3c dy 

d00„ 

7» v t — £. where c denotes the initial coordinate z of 
d-y 

the fluid particle, or the value of z whereby the vortic 
ity is considered equal to that corresponding to the mean 
motion at the same point. If a fluid element initially 
parallel to z keeps its orientation constant, then obvi- 

or is always constant, and therefore 

oc 



v i 3-?- = 0; but if the instantaneous velocity also has a 

3c • . 

component along z, 1 1 ~ 1 ' which represents the defor- 
mation of the element defined above along z itself, it 

may assume any values whatever and. therefore v' f-^- may 

o c 

be different from zero. Taylor shows that if I is suffi 
ciently small, in the special case of flow about a flat 
plate, the application of the principle of the transport 
of vorticity allows an equation to be obtained that is for 
mally identical with (5), derived from the entirely differ 
ent concept of Prandtl. 

5. Taylor's conception is certainly ingenious and 
may perhaps give useful results in those more complicated 
problems where the other theorems prove powerless. Never- 
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theless, it does not appear to the writer that it is neces- 
sary to have recourse to. a three-dimensional theory of tur- 
bulence in order to explain the limits of appl icability of 
formulas (7) and (7'), to give reasons for the erroneous 
results to which they may lead, and to substitute for 
these other formulas of greater generality. 

It is interesting to compare the results of the pre- 
ceding theories with those that may "be deduced by the cor- 
rect application of the equations of Navier, which are nat- 
urally assumed valid not only for the instantaneous motion 
but also for the mean motion. Under these conditions, 
there is obtained from the equation of projection of the 
momentum in the x direction 

p Qx dx By 

in which u and v are the instantaneous values of the 

velocity, or u = U + u' , v = V + v' ; u and _^v always 

being connected by the equation of continuity + = 0. 

dx dy 

There immediately results in place of (3) for the mean mo- 
t ion 



p ox dy dx 3x 



v m. 

By 



But if /u' 2 = I 9JJ it is possible to put u' = - V 

By By 



i — 



= a r 



3U 

By 

pendent of y, 



in which a may vary in time but is inde- 



v i 2fil - _ a A. 
3y 2 dy 



s ay 

du 



3U 




By 





and therefore, including as before the single constant a 
in 1 1 , 



_ * BP = _ 1 _EL 
p dx 2 3y 



3i , i i 

' 2 dx 



,s /au\ 

1 



If the motion is independent of x, (ll) reduces to 



_ I a_p = > i jl r I s 
P ax 2 dy [ W; J 



(11) 



(12) 
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which is formally identical with (7) and therefore, under- 
equal assumptions on theforra of I, leads to the same 
distribution of velocity. In particular, for the flow 
about a flat plate, 

dp 

and there is therefore obtained 

|S = u = k log y + b 

^y y 

On the other hand, the temperature distribution law 
for the same case is given by 

»il = ° = £(^)=£ O'gS) <»> 

which, by putting 



T = b U 



comes out identical with (12) (since 4™ is zero), and 

therefore assures the similitude of the diagrams of the 
velocity and temperature distribution in the turbulent 
flow about a flat heated plate. 

In the wake of a heated obstacle, however, at a great 
distance x downstream, equation (ll), with the same de- 
gree of approximation by which (9) is written, the torm 

^'jV J being negligible compared with |V 



gives : 



0 3x 2 3y 



V 



3 /au^ 2 



(14) 



which agrees with (9') resulting from the theory of Taylor, 
while the temperature distribution is always given by (10). 

6. From what we have said above, it therefore turns 
out that the lack of agreement of the experimental results 
with those deduced from (7) in the problem of turbulent 
diffusion in the wake of a cylindrical obstacle, is not a 
consequence of (4) and (6) and therefore of a possible in- 
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fluence of the local pressure gradient which, "by varying 
the quantity of motion of the particles , increases the co- 
efficient of heat transport with respect to that of impulse, 
as deduced hy Taylor, hut of the fact that the equation of- 
motion is not (9) hut (14). 

This conclusion does not, however, yet justify us in 
deducing any principle of general character as regards the 
possibility of the application of the Taylor and Prandtl 
concepts to the problems of turhulence. In fact, the fun- 
damental relations (4) and (6) whose validity we have just 
shown for the cases of the two-dimensional motion consid- 
ered, and which in the theory of Prandtl define the con- 
stancy of the momentum in the transport of turbulence, 
have "been deduced hy an entirely distinct procedure hy 
Karman, and by this method acquire an essentially differ- 
ent significance, which determines also the limits of ap- 
plicability. 

In 1930 Karman had already determined equations (4) 
and (6), assuming the condition that the disturbed motion 
corresponding to the turbulent agitation of the fluid par- 
ticles is statistically similar at all points of the field, 
differing from point to point only by the scales of time 
and of length. This assumption, as Karman observes, per- 
fectly corresponds to that which is normally made in the 
kinetic theory of gases and which permits the stresses due 
to the thermal agitation of the molecules to be simply ex- 
pressed hy means of local derivatives of the general ve- 
locity of the motion, and of the mean, molecular trajectory, 
and yields as a consequence the constant correlation between 
the components of the turbulent oscillation velocity, to 
which we have already referred ahove and which has been 
confirmed hy the experiments of Wattendorf (reference 11) 
and Kuethe at Pasadena, and by Eoichardt (reference 12) at 
G-ottingen (fig. 2). In this way, under the assumption 
that the general motion takes place very approximately 
along flow lines parallel to the X axis, there result the 
fundamental relations of Karman,-. 



which, being formally ident ical . with those estahlished by 
Prandtl, differ substantially hy the concept from which 
they were derived and by the definition of the length I, 
which in the theory of. Karman, is given by 
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I = 



dU/dy 

ii 

dy 2 



The meaning of (4) and (6) has been recently general- 
ized "by Dede'bant , Schereschewski , and Wehrle (reference 
13) , who have shown that the Karman relations may he ob- 
tained by applying the law of similitude not to the con<- 
figuration of the velocity field, but only to the law of 
distribution of the disturbed velocity; that is, by assum- 
ing that the disturbed velocity distribution functions be- 
come identical by a suitable change of scale. 

It should still be observed that the equation (3') 
leads to the same laws as those derived from (12) by assum- 
ing that the fluctuations of the vorticity at each point 

are proportional not to =^ but to ^~ I ^~ , and 

therefore not to the derivative of the mean rotation but 
to the derivative of the displacement, 

7. Tor the more accurate comprehension of the effec- 
tive meaning of the preceding formulas and of the real pos- 
sibility that the concepts explained above offer for a sat- 
isfactory solution of the problem of turbulence, it is of 
great aid to consider the motion along flow lines of non- 
zero curvature. For simplicity of treatment and by analo- 
gy with what has been done above, we shall suppose that 
the mean flow always takes place in a plane and along arcs 
of concentric curves. We shall denote the velocity tan- 
gent to the circle of radius r by v^. and the radial ve- 
locity by v r . The natural extension of the concept of 
Prandtl (references 14, 15, and 16), leads to the conclu- 
sion that the fluid particles in being transferred from 
layer to layer as a consequence of the turbulent agitation, 
maintain their velocity moment constant with respect to 
the center of rotation - that is, T{ r = constant = c. 
A fluid particle, therefore, which arrives at a layer of 
radius r after transversing a radial distance I has 

c t -V t t 

decreased its own velocity by Av + = - — 5- I = — ~ I , 

r dv 

and since it had initially an excess of velocity l f — , 

dr 

the increase in velocity that it has with respect to that 
of the layer to which it arrives, becomes 
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v' 



(15) 



The mean increment in the momentum communicated to the el- 
ement of radial thickness dr and of length r d cp, is 
therefore, 



d"M = r d<p dr v« r l» 



_d_ 
dr 



(v t r) 



to. which there corresponds a virtual tangential stress 

dr (V * r) 



T = p 



V' V d 



(16) 



different from zero on account 
in general exists between v' r 
and l'. 



of the correlation which 



and v' t , or "between 



v 1 



There is thus obtained the equation of motion in the 
mean direction of flow under the assumption that v^ is 
constant along a flow line of the mean motion: 



lid' 

( T 

p r dr 



p dcp r dr L 



v« r l' 



d f 
d7 (Vt 



r)3 (17) 



Prom (17) and (16), however, there follow erroneous 
results or results that have not been confirmed by experi- 
ments. Karman, in fact, observes that from (16) it would 
be deduced that T becomes zero at the point at which 
v-fcr assumes its maximum value, and that the sign of T is 
determined by that of ~ (v^r) . Now Wattendorf (refer- 
ence 17) has recently conducted tests on a channel with 
circular axis, having an elongated straight rectangular 
cross section (ratio of the sides l/l8) , so that the ef- 
fect of the secondary flows might be assumed negligible, 
and from the fall of pressure along the channel, and from 
the direct measurement of the stress tangential to the 
walls, deduced the diagram for T along the radius. Karman 
remarks that from this the relations given above between 
T and v t r (fig. 3), are not confirmed. However, this 
observation of Karman and Wattendorf does not appear suf- 
ficient to the writer to invalidate the assumption of 
Prandtl, and therefore, indirectly,' the hypothesis of con- 
stant moment of velocity during the turbulent transport. 
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In fact, from the equations of Navier, applied to the 
particular type of flow now considered, there is obtained 
for the instantaneous motion: 

1 BP Bv t , /v r l 3v t \ 

- — = V r - + Vf ( — + — J 

p r3<p Br * V r r 9?/ 

and for the mean motion 



which, with the assumption of Prandtl, becomes 

Now since for the determination of T, Wattendorf availed 
himself of the relation 

I ± (r a T) = 3P 
r dr V ' Qcp 

By comparing (19) with (17) , it follow that the virtual 
stress which gives rise to the pressure gradient is not 
that which results from (16) but that given by 



T = — -~ 
r 3 



/ r v' r ~- [l* v t r)] dr + a 



A very remarkable observation on the theory of Prandtl 
is also made by Taylor, who remarks that the theory of the 
transport of the moment of momentum necessarily leads to 
expression (16) for the apparent stress, and therefore a 
rotation of the entire system as a rigid body with angular 
velocity Q would change the tangential stresses, increas- 
ing them by an amount 2p v' r I ' Q. 

Taylor's theory does not present this incongruence 
inasmuch as the addition of a constant vorticity to the 
entire fluid field does not have any effect on the trans- 
port of vorticity itself. 



8. In this respect, too, however, it seems fitting 
to show how every difficulty for the determination of the 
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expression for the fluctuation velocity might he overcome 
not hy assuming a certain transport length, and more or 
less plausible hypothesis based on the law of motion of 
the particles during their transfer from layer to layer, 
but by applying the principle of statistical similitude 
according to the concept of Zarman mentioned above. 

In the coordinates r, <p of the plane , the equation 
of the transport of the vortices becomes: 

^ D^ + v t |- Dt + r v r D\|r = 0 (20) 

where t denotes the flow function of the field and D the 
symbo 1 

3 r 2 r Br r 2 3<?> 2 

Setting ty x the flow function of the motion corresponding 

to the turbulent agitation, and denoting by Vt and <d , 
respectively, the velocity and the mean vorticity at any 
point, equation (20) may be written: 



Let us now make the assumption that with respect to a 
system of axes with origin at any point P of the field 
and moving with P with the same velocity, the field of 
flow in the neighborhood of P may be considered as sta- 
tionary. We can then put, in the neighborhood of P 



0) 



*Vt\ 
dr L 



▼t p 

(r - r p ) --~ (r - r p ) + ... 



-» P+ (g) . (r-rp) + ... 



> (22) 



since we should evidently put in the second member of the 
first of equations (22) the disturbing velocity of the 
fluid element. 



There results immediately if-, in analogy to what was 
already done by Karman for the rectilinear motion, we put 
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r - r p = I r x ; ti(r,<P) = Af(r x ,cp) 

9t: 



rp Vdr /p 



(r ~ r p ) g^Dt, - 3 ~- + 



Br 3<p rx Bqp 8r 



and 



Vdr/p r p 



8 2 /doa\ 3f 



i ar x a<p ^ a<p 3r x 



(23) 



If the form of the function f which defines the 
fluctuations of the velocities due to the turbulent motion 
is to be independent of the particular position of the 
point P, then it is necessary that 



dr /_ r. 



>dV 
V 



P . 



/d(0\ 



A 



Vdr/L ~ 1 



or the characteristic length I becomes 



I = 



dVt 
dr 



It 
r 



d_U> 
dr 



(24) 



while A is given by 

A = I 2 



/dV 



_Lt\ 
dr J 



It] 
r 



(25) 



and therefore: the amplitude of the turbulence velocity 

V 

is proportional to I 



Vdr 



- ) whereas Prandtl put the 
F_Z 



amplitude of the fluctuations themse^ves_projpo_rt_ional to 

Vdr r J 
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Since the first expression in a rigid, rotation is zero, it 
has no effect on the amplitudes of the oscillations them- 
selves, and therefore on the expression for the virtual 
tangential stress. 

9. It is interesting to compare (24) and (25) ob- 
tained by the principle of similitude of Karman from the 
equation of the transport of vortices with those that may 
he deduced by extending the procedure already mentioned of 
Dedebant, Schereschewski, and Wehrle to the type of flow 
we are now examining. 

Let us assume that the distribution functions of the 
perturbed velocities at each point of the field may he 
made to become idential simply by a suitable change in the 
scale of the velocity and of the time, and let us suppose, 
following the method indicated by De'debant , Schereschewski, 
and Wehrle, and moreover, the procedure used by lorentz 
(reference 18), and by Chapman (reference 19) in the kinet- 
ic theory of gases, that the distribution function f is 
very nearly the same as that of Maxwell, so that indicat- 
ing the latter function by f 0 , (f Q = a e"*•' b ^ n, +v ' >), 
f = f 0 (1 + €), in which £ is a small quantity of the 
first order. Assuming as unit velocity at each point of 
the field the mean quadratic variation a of the disturbed 
velocity, the values of f are everywhere identical; the 
values of € should therefore be the same. Now e, which 
it should be possible to represont by a series in the com- 
ponents of the disturbed velocities u 1 and v' (which 
we shall now assume as referred to the mean value a as- 
sumed as the unit) , contains the terms which define the 
nonunif ormity of the field; or those of the velocity of the 
general motion and those which De'de'bant calls the "densi- 
ties." In fact, in order to make the statistics of the 
disturbed velocities comparable, it is necessary that the 
time of observation of the velocities themselves vary from 
point to point in such a manner that the number of the 
fluctuations examined be everywhere equal. Now in time t 
the particles observed are proportional to 1? ; resulting 
in a density, if N is the constant number of fluctuations 
considered proportional to N/ 1 8 . The corresponding non- 
uniformity is therefore now defined by the variation of t. 
The distribution function, however, is an invariant; that 
is, it does not change in form or in value. Since f Q is 
by itself invariant, € should be the same; but the invari- 
ant quantities, functions of u' and v 1 and of the deriv- 
atives along the same axes of t and of the general veloc- 
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ities of the components V and V are functions of the 
following elementary invariants: 

u« + v' at 
3x By 

corresponding in polar coordinates and by the, assumed char- 
acteristics of the motion considered to v' r f-t, and 

dr 

u ,» au + v ,« ai + u , v , (as + an 
ax ay vay ax / 



corresponding to 



, , /ATt V t\ 



There results, according to Chapman, the following ex- 
pression for €: 



/dV t V t \ dt 



in which A and B, on account of the homogeneity of the 
formula, ^'t' v, r being simply numerical, should be put 
proportional to t and a, respectively. Since e should 
be independent of the coordinates of the point at which 

the distribution function is determined a — — = constant 

dr 

or a is proportional to and t f— — '. - -— ^ = con_ 

1 r dVt 
stant , or t is proportional to -T™ — • Bnt r~~ " 

dYj V t dr 

dr r 

V+- 

-- is equal to 7 ,„ and therefore, the mean quadratic 
r r< r 

variation of the velocity of fluctuation is 



where 1 is now given by 
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There is thus found the property that the amplitude 
of the velocity fluctuations due t o the turbulent agita- 
tion, is proportional to y r( ^ "but the characteristic 

length I does not depend, as before, on the derivative 
of the rotation, hut still on the derivative of the 7 r qj« 

This conclusion immediately makes the results of the pre- 
ceding investigation more understandable: When the rota- 
tion is zero at every point, as in the velocity field 
which is generated in a viscous fluid about a rotating 
cylinder of infinite axial length, it is logical to assume 
that turbulence cannot take place. 

Now while (26) and (27) do not tell us anything about 
this point, (24) and (25) assure us of the impossibility 

of a motion having the assumed characteristics since — ~ - 

dr 

0 and therefore I = <». 



10. The considerations just developed bring out 
the difficulty of developing according to the methods of 
classical mechanics any theory whatever based on the con- 
cept of transport, inasmuch as the conclusions to which 
they lead are intelligible only insofar as the relations, 
which the introduction of the concept of mean distance I 
permit to be written down, may be interpreted by means of 
the principle of similitude of the turbulent oscillations, 
while a generalization of the relations themselves with- 
out this check easily leads to erroneous conclusions. 
Thus in the case of plane motion with rectilinear stream- 
lines (4) and (6) are correct since dU/dy represents the 
excess of velocity with respect to the fluid layer at 
height y of the particle distance I, as well as the 
fluid displacement between the two contiguous layers, 

while (15) leads to unintelligible results since ^Zt + Y_L 

dr r 

is proportional to the excess of moment of momentum while 

dv-t vt 

the fluid displacement is given by j-— - — — . On the 

other hand, the reasons for such difficulties are easily 
understood, for although the concept of the mixing length 
is very intuitive it is not" prescribed for a fluid parti- 
cle which penetrates a certain layer after having trav- 
ersed more layers, how many of the dynamic and kinematic 
characteristics that it possessed initially, it conserves 
after each crossing, 

11. With the problem of the transport of momentum 
and of vorticity is intimately connected the theory - re- 
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markable in many respects - developed recently by G. D« 
Mattioli (reference 20). The latter supposes that in tur- 
bulent agitation each fluid particle maintains constant 
not only its momentum, as Prandtl. assumes , but also its 
rotation, as assumed "by Taylor. Mattioli, who considers 
the motion in tubes of circular cross section, derives the 
following equations: 

where V is the general velocity, e = I 2 k, in which I 
has the meaning' given above, while k is a function which 
has the dimensions of a frequency and is called by Matti- 
oli the "mixing frequency." 

How if Prandtl's theory is applied to this particular 
type of flow, we have, using the above notation 

L i?. r = ± (T^ ^) 

p dx dr V r dry 

which is identical with the first of equations (2.8) putting 



V v' = I 2 k (29) 

On the other hand, from the equation of Havier, which rep- 
resents the motion in the direction of the tube axis, we 
have: 

l p. = v^ r 
p dx 

where to* is the instantaneous fluctuation of the vortic- 
ity at any pdint." Therefore, according to the Taylor con- 
cept, 

= v 

dr 2 

from which 

p dx dr P dx dr 



or 



d 

dr 



L / d 2 V\ 1 dp 
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Equation (30) does hot give satisfactory result s ,- how- 
ever. Taylor, as explained above, concludes therefrom 
that the theory -of t ransport- -of- -vort-icity could not "be ap- 
plied in this case except by considering the turbulent agi- 
tation as it is actually, namely, three-dimensional. Mat- 
tioli, whose theory is one-dimensional, insofar as the 
transport of the masses takes place only with radial tur- 
bulent velocities, and the continuity of the mass itself 
is restored by associating with the discontinuous turbu- 
lent transport a continuous transport that is. always radi- 
al, does not at all consider (30) or the connections set 
up by the equations ofUavier, and since he shows that the 
first member of (30) represents the increment of vorticity 
communicated to an element of fluid, he deduces that, for 
the permanence of the motion, there should be applied to 
the same element a corresponding couple which he puts 

equal to ar *■ • 
dr 

He maintains (reference 21) that this couple is due 
to the viscosity of the fluid and therefore dissipates at 
each point the increment of vorticity which the turbulent 
agitation produces. This interpretation is open, however, 
to some reflection even if the fact that the fluid, in the 
original theory, was assumed to be perfect. 

The first member of the second of equations (28) , 
which I shall denote by A, represents the excess of vor- 
ticity which is to be dissipated in the time t 0 corre*- 
sponding to the mean period of the oscillations. Now the 
velocity of dissipation of the vorticity « in a fluid of 

kinematic viscosity u is given by ~ = uAu) , by 

— cL t — - — 

which the mean velocity f^. in time t n may be put f™ - 
_ dt . ° dt 

vA5 and the dissipated vorticity will be 

jjSL t 0 = V t 0 AW = A (31) 
at 

In this way it may be understood why in the second 
member of (28) there does not appear the viscosity, since 
putting the time t Q necessary for the dissipation of A 
proportional to l/u seems logical, and in any case agrees 
with what may be deduced in several simple cases (refer- 
ence 22). Equation (31), however, would make A depend 
not on the second derivative of V but on the third de- 
rivative. 
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Apart from the difficulty of justifying the second of 
equations (28) which has already been the subject of some 
remarks by the author (reference 23), the theory of Matti- 
oli leads to results which agree both with those obtained 
experimentally and with those deduced from the theory of 
Prandtl and Karman; in particular as regards .the logarith- 
mic law for the resistance of smooth tubes as indicated 
by Mattioli himself (reference 24) and the logarithmic dis- 
tribution law of the velocity at the wall of the obstacle, 
as may easily be deduced. 

11. The problem of the transport of momentum or vor- 
ticity in turbulent flow, presents a particularly sugges- 
tive aspect and is susceptible of a general solution if 
the phenomenon is considered as belonging to the domain of 
statistical physics and is therefore studied with tho meth- 
ods appropriate to the latter. In this connection, it 
should be observed that the statistical character of the 
flow is assured, not only by the irregularity of the ve- 
locity fluctuations, but by the existence of universal 
laws of distribution of velocity independent of any ini- 
tial condition of the flow. Important progress in this 
sense is represented by the theory recently developed by 
H. G-ebelein (reference 25). 

G-ebelein assumes that in the fluid motion, the circu- 
lation of the momentum, and the energy diffuse in space 
according to the same law of probability that governs the 
diffusion of matter in statistical phenomena; that is, in 
which the motion of any particle from point to point is 
not determined in an unequivocal manner, with certainty 
(as in problems of deterministic mechanics), but takes 
place according to a law of probability, which is a time 
and space function. Under these conditions, the density 
of the fluid in a space element, which is proportional to 
the probability with which molecules of fluid are found 
in the element considered, received in each elementary 
volume an increment p div. V, which by the continu- 

dt 

ity of motion should be zero, and whose analytical ex- 
pression is given by the equation of Kolmogoroff: 

k=3 

i=3 a 
2 (b ik p)=0 (32) 

i=l 9yi3y k 
k=l 



dp 

it + p 



div. V 



8t 



i=3 
2 



2 -S-(uip). 
i=l 3yi 
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where y^ now represents. any of the three coordinate axes 

(y x , y 3 , y 3 ) t ni the component in the y^ direction of 

the mean velocity of the motion, while hit depends on 

the mean quadratic variations of the velocity itself, 

Ho re precisely, let P (x lf x 2 , x 3 , y a , y 2 , y 3 ; t; t .+ A ) 

he the function which defines the probability with which a 
particle, which at time t is at the point whose coordi- 
nates are x t , x 2 , x 3 will he at the point & , y 2 , y 3 at 
time t + A, and Ai and Bit, the moments of first and 
second order, respectively, of P, or the mean values and 
the quadratic variations, respectively, of the distribu- 
tion function P, that is 

+CO +00 +00 

A i = / / / (yi - x i) p d yi d y 3 d y 3 

—CD -00 —00 
+00 +03 +00 

B ik = / / / (yi - x i) (yk - x k) p d yi d y 3 d y 3 

—00 — 03 —OS 

Then u* = lim ; b ik = lira 

1 A=0 * A=0 

Thus the equations of motion, according to the statis- 
tical mechanics, do not require a knowledge of the form of 
the function P, but only of the values of the static mo- 
ments aid of the second order moments of P itself, pre- 
cisely as in the deterministic mechanics, in which the gov- 
erning equations of motion contain only the static moments 
of inertia of the mass in motion and do not require a knowl 
edge of the form of the distribution law of the mass. 

According to the hypothesis given above by Gebeloin, 
if the mechanical characteristics are diffused by the 
same law as the density, the increment of the component of 
the momentum of an element . dT in the direction yi , 

which, on account of the translation equilibrium along the 
same axis, is equal to - d^/By^ is given by 
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-at" + 2 ; §57 " 2 ay-ay^ OikP^i) = - dy - 

while analogous equations are obtained for the directions 
y 3 , y 3 .. 

If, instead of the momentum, it is assumed that it is 
the vorticity that is diffused, there is obtained: 



Byi 



Bu^ 

pU2 557 



6uj_ 

Sy 2 ~ 



12. The solution of the problem requires a knowledge 
of the coefficients b-^ and this is obtained by G-ebelein 

by using a relation between the mean quadratic variation 
of the velocity and the characteristics of the mean mo- 
tion, and the expression for the "characteristic time" 
T Q ( Verweilzeit) which is a measure, so to speak, of the 
statistic or deterministic character of the phenomenon un- 
der consideration, in the sense that the phenomenon itself 
observed for an interval of time less than T 0 appears to 
be governed by a deterministic law while observed at inter- 



vals of time very large compared to 
tistic phenomenon. 



appears as a sta- 



G-ebelein' assumes as a fundamental theorem, that the 
mean quadratic variation of the perturbed velocity is pro- 
portional to the fourth root of the vorticity, correspond- 
ing to the mean motion 



4 i 

a = yco 



(33) 



but this proposition is, in fact, not demonstrated, at 
least not in a convincing manner, and therefore appears 
essentially as a hypothesis to which G-ebelein gives an ex- 
perimental confirmation based on the tests of Hikuradse , 
of 1926. Since, however, the determination of the mean 
quadratic variation in these experiments is not made by 
direct measurement, it seemed proper to the author to ver- 
ify the hypothesis of Gebelein by means of the results 
obtained experimentally with a hot-wire anemometer by 
Reichardt and by Wattendorf in a very elongated rectangular 
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channel so as to approximate as far as possible two-dimen- 
sional motion. The results which have "been calculated 
using the values indicated above, are shown in figure 4 in 
which is also indicated "by dotted line the theoretical 
diagram according to the assumption of G-ebelein. The ap- 
proximate character of (33) is deduced "by the comparison. 
On the other hand, if the motion takes place with a varia- 
ble pressure gradient in the direction of the motion it- 
self "between a convergent and a divergent section, (33) 
does not give any indication of the eventual dissymmetry 
between the mean quadratic variations of the components 
u» and v' of V, while it is known, for example, that 

in the motion within a convergent yy/u 1 2 decreases, along 
the axis of the channel, abo\xt in proportion to L, with 

which the mean velocity increases, while J v 1 2 increases 
in proportion to L 1 2 . In this connection, it is very 
desirable that systematic tests be conducted to determine 



the law of variation of the dependence of J u ,2 ~ and of 



J 



v' on the mean motion, the knowledge of which is es- 



sential for any theory of turbulence. 

The characteristic time T 0 is determined by Gebelein 
after an actual analysis of the reasons which could give 
rise to a statistic phenomenon destroying the causes which 
tend to produce it in a deterministic manner. The cause is 
essentially the same as that leading to the production of 
vortices at the contact of rigid walls which, either on 
account of the roughness of the walls themselves or the 
disturbances which may be produced, for example, at the 
mouth of the tubes, takes place not according to a deter- 
minate law but by pure chance. This is what makes the 
vortices in the boundary layer, in two-dimensional motion, 
have only in the mean a direction parallel to the walls 
and normal to the mean motion; to this mean vorticose 
layer is added a layer in which the axes of the vorticese 
are disposed along any direction whatever. A comparison 
of the results to which the theory leads with those ob- 
tained experimentally would indicate that the axes of the 
perturbed vorticese are essentially directed along the 
mean motion, and this conclusion appears in singular agree- 
ment with the hypothesis of Prandtl already referred, to, 
namely, that the. transport of particles in bound turbu- 
lence takes place- by just the vortices having . this dispo- 
sition. 
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13. The calculation of T 0 has "been made by Gebelein 
for several- simple cases (flow about a flat plate, within 
a convergent and a divergent tube, and within a straight 
tube)-; and some of the results obtained by him by relative- 
ly simple methods, although at the price of over-simplified 
assumptions, we shall now point out in what follows. 

In the problem of the flow about a flat plate, the 
equation of Gebelein obtained by assuming that the momen- 
tum is diffused according to a law of probability corre- 
sponding to that of the diffusion of the mass, considering 
the fluid as perfect, becomes 



1 dT a 2 



( v ' 2 T o u i) = 0 (34) 



p d y 2 d ^ 

while the diffusion of the vorticity leads to the equation 



1 dT d / a du x \ 

~ T— = J — (v 2 > T 0 - — — ) = 0 (35) 
P dy 3 dy 2 V 0 dy 2 / 

2~ 

By the assumption made on the relation between v' 
du, 

and — — and with the expression for T 0 calculated by 
Gebelein, we have for the two cases, respectively: 



d 2 / duj, \ 

d"?7 7% I 7 7 Ul) = ° (36> 



_d_ 



= 0 (37) 



Equation (36) gives u (y g ) = J a + b log y~ , while 
from (37) is obtained 

u i(y 2 ) = a + b log y 2 (38) 

which is the well-known logarithmic formula confirmed by 
the tests of ETikuradse, mentioned above. 

Thus, it is the vorticity and not the momentum which 
is diffused according to the same law of probability as 
that corresponding to the diffusion phenomenon of the mass. 
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This conclusion may, however, be invalidated by the 
assumptions which are at ' the base, of the . calculations; in 
particular, it may be difficult to admit that the mean 
quadratic variations of the fluctuation velocities for the 



flow about a flat plate are proportional to 

1/4 



du x 

ay" 



therefore by (38) to y" A/ "*. This even appears to be con- 
trary to the principle of the constant correlation between 
u' and v', which would lead,, in this problem, and in 
accordance with the principle of Karman, to the equation 



v' = constant 



(39) 



It may, nevertheless, easily be seen that the con- 
clusion stated above holds true even if equations (39) 
are admitted. In fact, analogous to what is done in the 
kinet ic . theory of gases, it is possible to assume 

I 



if I has the meaning of "the mean free path tt already 
considered many times above; there is deduced 



_d_ 
dy 




o r 



- 0 



Sku 
dy 



= constant 



which leads to (38) , as already previously derived, 

14. The theory of Gebelein, developed according to 
the methods of statistical mechanics, thus allows the af- 
firmation of the principle of diffusion of the vorticity 
in turbulent .motion. This principle had been confirmed by 
Taylor, but its application, using the methods of classical 
mechanics, had not led to, a sat isf acto^y solution.. The 
real reason for this fact appears to depend On the circum- 
stance that a theory based on the concept of transport 
and developed according to classical mechanics, would re- 
quire 'a knowledge of the history of the particle during 
the transport, whereas statistical mechanics, observing 



28 



N.A.C.A:. Technical Memorandum No. 799 



the phenomenon at intervals of time for which every effect 
of cause is destroyed, renders itself independent of any 
knowledge of the states assumed "by the particle during its 
deterministic motion. 



Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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(a) Wattendorf (C.I.T.) 

(b) Reichardt (Gottingen) 
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Figs. 3,4 
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